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Topological or deconfined phases are characterized by emergent, weakly fluctuating, gauge fields. In con- 
densed matter settings they inevitably come coupled to excitations that carry the corresponding gauge charges 
which invalidate the standard diagnostic of deconfinement — the Wilson loop. Inspired by a mapping between 
symmetric sponges and the deconfined phase of the Z2 gauge theory, we construct a diagnostic for deconfine- 
ment that has the interpretation of a line tension. One operator version of this diagnostic turns out to be the 
Fredenhagen-Marcu order parameter known to lattice gauge theorists and we show that a different version is 
best suited to condensed matter systems. We discuss generalizations of the diagnostic, use it to establish the 
existence of finite temperature topological phases in d > 3 dimensions and show that multiplets of the diagnostic 
are useful in settings with multiple phases such as 1/(1) gauge theories with charge q matter. [Additionally we 
present an exact reduction of the partition function of the toric code in general dimensions to a well studied 
problem.] 



I. INTRODUCTION 

There is currently much interest in condensed matter sys- 
tems that exhibit ordering captured by a gauge theory lack- 
ing a local order parameter. Such systems offer a contrast to 
the classical broken symmetry paradigm built around the no- 
tion of an order parameter introduced by Landau and these 
phases said to be topological or deconfinecU*£ as they also ex- 
hibit quasiparticles with fractional quantum numbers 2 . These 
phases are of interest in the study of strongly correlated quan- 
tum systems — where the quantum Hall states- and resonating 
valence bond liquids^ are the canonical examples. Their in- 
terest has been further enhanced by Kitaev 6 and Freedman'si 
proposal of utilizing their subset possessed of non-abelian 
braiding statistics for quasiparticles, for the construction of 
a quantum computer. 

The theoretical description of such phases is in terms of a 
gauge theory, possibly with low energy matter. In the simplest 
cases the low-energy theory is a purely topological gauge the- 
ory such as Chern-Simons theory in the case of the quantum 
Hall states- or the BF theory in the case of resonating valence 
bond liquids 8 -. Indeed, the term topological phase dates from 
these early instances; absent a better standard term, we will 
use it here to refer to all phases with an emergent gauge field. 
The explosion of recent interest has come from the construc- 
tion of a wide variety of lattice models that realize a variety 
of long wavelength theories including those with dynamical 
gauge fields^. 

In idealized models the existence of such phases is transpar- 
ent for one can readily show that the gauge field: (a) exhibits 
weak fluctuations and thus a perimeter law for the Wilson loop 
in ground states, (b) mediates a non-confining force between 
its sources ("quarks"), and (c) [where the low energy theory 
is purely topological] leads to a topology dependent ground 
state degeneracy. For more realistic Hamiltonians, unfortu- 
nately, this transparency is lost. The problem with the first 
two characterizations on our list is that, as is well known from 
the lattice gauge theory literature, they fail in the presence of 
dynamical sources ("matter") for the gauge field. 



In the condensed matter setting, the gauge field is only 
emergent and necessarily comes with higher energy excita- 
tions that carry a gauge charge, making this a generically fa- 
tal flaw. The situation only gets worse at finite temperatures 
where such matter excitations must enter the statistical sum. 
The third characterization can continue to work in strictly 
topological phases precisely at T = but lacks a clear mean- 
ing beyond that limit and outside that subset of deconfined 
phases^ 

All of this is an unsatisfactory state of affairs. As in Landau 
theory, it would be nice to have a fixed time diagnostic oper- 
ator that one can compute when handed a ground state or a 
Gibbs state to decide whether it exhibits the ordering charac- 
teristic of a given topological phased 1 In this paper we report 
the construction of such a diagnostic which teases out the un- 
derlying weakly fluctuating behavior of the gauge field. We 
do so via a detour into the theory of sponge phases where our 
diagnostic has the interpretation of a line tension — picking 
up a line of analysis begun by Huse and Leibler= a while 
back. Remarkably, this diagnostic turns out to be a space- 
time generalization of the so-called Fredenhagen-Marcu order 
paramete r 13 i 14 known to lattice gauge theorists. 

In this paper we also begin the process of applying these 
ideas to a wide class of systems. While the bulk of our pa- 
per is concerned with the Z2 gauge theory with matter (or Ki- 
taev's toric code 6 with generic perturbations) we sketch the 
generalization to different low energy gauge structures includ- 
ing cases where it is necessary to use a multiplet of diagnos- 
tics. Notably, we use the diagnostic to establish the survival 
of topological phases to finite temperatures in d > 3-^. 

Turning to the contents of the paper, we begin in Sections 
II and III by reviewing the formulation of the Z2 gauge theory 
with Z2 matter and its reformulation as a statistical mechan- 
ics of surfaces with edges. This leads us to the identifica- 
tion of the line tension diagnostic in Section IV and its oper- 
ator formulations in Section V. In Section VI we discuss why 
condensed matter systems prefer a particular operator formu- 
lation. Sections VII and VIII deal with applications of the 
ideas to finite temperature phases and more complex phase 
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diagrams. The conclusion flags some open questions. Finally, 
two appendices establish some useful results on Kitaev's toric 
code: firstly, we obtain its partition function in general di- 
mension via a reduction to a classical pure Z2 gauge theory; 
secondly, we demonstrate perturbative stability of its topolog- 
ical orde r 16 ! 17 , and find it to persist to finite temperature only 
for d > 3, in accordance with Ref. [18]. 

Before proceeding, we should point out an important an- 
tecedent to our work which also addressed the challenge of 
exhibiting topological order away from the idealized mod- 
els. Hastings and WeniS gave a T = continuity construc- 
tion starting with idealized Hamiltonians that generates a 
Hamiltonian-dependent gauge field operator which exhibits 
a perfect perimeter law (with zero coefficient) for the Wilson 
loop in deconfined phases. By contrast, our construction is 
Hamiltonian-independent in the spirit of an order parameter. 



n. Z 2 GAUGE THEORY WITH MATTER 

Consider the Hamiltonian of the Z2 lattice gauge theory, 

p ledp l I sedl s 

supplemented by the constraint that we restrict its action to 
"gauge invariant" states defined by 
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where the gauge (oj) and matter (t' v ) operators act in spin 1/2 
Hilbert spaces that live on the links / and sites s, respectively, 
of a (hyper) cubic lattice in d dimensions. The subscripts s, 
I and p denote sites, links and plaquettes. dp and dl are the 
boundaries of the corresponding objects. 

For r = / = our Hamiltonian is equivalent to Kitaev's 
toric code as noted in the original paper 6 In that case the de- 
generate 2 N " C ground states (N nc is the number of independent 
non-contractible loops on the lattice) are easily constructed 
and exhibit a perfect "perimeter" law for the contractible Wil- 
son loop 



W[C] = (Y\o) 
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Consistent with this ground state description, one sees that the 
excited states contain non-interacting charges, free to sit at ar- 
bitrary locations in this particular model. The remaining part 
of the excitation spectrum is described by purely gauge exci- 
tations. These are vortices of the gauge field, or visons, which 
are isolated plaquettes for which W [C] = — 1 . Both charged 
and gauge excitations are separated from the set of ground 
states by a finite gap. 

The full T = phase diagram of the Hamiltonian ([]]) has 
two distinct phases. The first is the deconfined phase that ex- 
ists for small perturbations of the Kitaev point. (For a more 
detailed analysis of the toric code phase diagram, including 
T > 0, see App. [A]) The remaining phase is the confined- 
Higgs phase which takes its name from the limiting regions 



that were shown to be smoothly connected by Fradkin and 
Shenkeri£. 

How do we distinguish between the two phases? Unfor- 
tunately, except along the line J/Tm = 0, the Wilson loop 
W [C] does not readily detect the phase transition — it exhibits 
a perimeter law everywhere (and also at T > 0). As previously 
remarked, the potential between two test charges is also not a 
sharp diagnostic due to the presence of dynamical matter. 



III. SURFACES, EDGES AND THE SYMMETRIC SPONGE 

In order to get to our diagnostic, let us now reformulate the 
Hamiltonian gauge theory considered above as the classical 
statistical mechanics of a system of membranes. To do this 
we first write down the path integral formulation with a dis- 
cretized imaginary time. This is governed by the action^: 
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where the O/ are the gauge fields living on the links, and T s 
the matter fields residing on the sites of a d + 1 dimensional 
hypercubic lattice, with d > 2 so that we can discuss a de- 
confined phase. In writing this spacetime symmetric form, we 
have moved away from the anisotropic limit that is needed for 
precise reproduction of the Hamiltonian problem but which is 
not important for our purposes. The reader should also bear in 
mind that K and J refer to different quantities, with different 
units, in the Hamiltonian (fl} and the action ©. However, as 
their physical import is very similar it is useful to keep this 
notation. 

We next rewrite the partition function as: 



Z (K,J) = (coshK) N » (cosh/) 
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where N„ / are the number of plaquettes and links, respec- 
tively, in the lattice. Multiplying out the products in the brack- 
ets and performing the trace annihilates all terms in which any 
Tj,G/ appears an odd number of times. For 7 = 0, this allows 
only for the presence of closed surfaces, a surface being de- 
fined as containing a given plaquette p if the factor of ri/e3p ®i 
appears in the sum. Switching on a non-zero J allows the sur- 
faces to have free edges. Thus we have rewritten our gauge 
theory as a statistical mechanics of surfaces with edges as 
promised. In this formulation each area element of the sur- 
face has a bare surface tension tanh/C and each edge a bare 
line tension tanh/. The actual geometrical properties of the 
different regimes and phases are determined by the renormal- 
ization of these quantities due to the entropy of the surfaces 
and edges. 

The phase diagram of the problem is sketched in Fig. Q] and 
includes more structure than is significant for the gauge theory 
as we explain below. We direct the reader's attention to the 
following features. 
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FIG. 1: Schematic phase diagram of the Z2 gauge theory coupled to Z2 matter (4), interpreted in analogy to sponge phases (see text). The 
inscribed cartoons capture the behavior of the Huse-Leibler horseshoe (|6j and allow detection of an effective line tension. 



Pure gauge theory: As already noted, for / = 0, the system 
excludes edges. All surfaces are closed, which allows a con- 
sistent definition of an inside and an outside. For small K, the 
total surface area remains small, and the volume of the inside 
of the surfaces is tiny, whereas the outside fills most of space. 
In this regime there is a macroscopic surface tension which is 
in one-to-one correspondence with an area law for the Wilson 
loop, as the reader can check by tracking the latter through our 
mapping. 

At a larger K = K g , the inside of the surfaces percolates 
as well. This, first, transition is purely geometrical and thus 
does not produce thermodynamic or local operator singular- 
ities, e.g. the surface tension is analytic across the bound- 
ary. In the language of membranes, the percolating phase is 
called the asymmetric sponge, as inside an outside have un- 
equal volumes. The transition is dual to the one encountered 
in the Ising models in d = 3, where the minority phase starts 
percolating at a temperature below the thermodynamic T c . 

At K c > K„, there is a thermodynamic transition which cor- 
responds to deconfinement in the pure gauge theory. This is 
signaled by the vanishing of the macroscopic surface tension: 
In the Wilson loop acquires a perimeter law, while for the 
sponge the inside and outside volumes become equal whence 
the new, deconfined, phase is labeled a "symmetric sponge". 

Regime with heavy dynamical matter) 3 - 2 - Turning on / 
causes edges to appear. Even at small J, where the bare line 
tension is large, edges have the qualitatively important effect 
of causing the macroscopic surface tension to vanish whence 
the Wilson loop exhibits a perimeter law for any K and / > 0. 
However, the phase transition between the asymmetric and 
symmetric sponges is not destroyed at small J even if the la- 



beling is now problematic as inside and outside are not dis- 
tinct. There continues to be an intuitive distinction in the ten- 
sion associated with surfaces between the two phases as illus- 
trated in the bottom panels in Fig. [2] In the symmetric sponge 
phase, the Wilson loop on the right follows a perimeter law 
despite the large area present, whereas on the left, a perimeter 
law arises because dynamical matter allows the surface area 
to be proportional to the perimeter of the loop by introducing 
an edge that runs along it. 

Regime with light dynamical matter: As the bare line ten- 
sion, tanh7, is increased further, edges proliferate. Of most 
interest to us, this has the effect of terminating the symmetric 
sponge phase and driving the system into the "sponge with 
free edges" which is the Higgs region of the gauge theory 
in the present formulation. This transition is most vivid at 
K = °°, where surfaces drop out and we are left with the sta- 
tistical mechanics of loops alone. In fact these loops the high- 
temperature expansion of the Ising model, with a vanishing 
line tension at J c signalling the onset of ferromagnetic order, 
a thermodynamic transition which separates the deconfined 
from the Higgs phased 

The sponge with free edges is accessible, without crossing 
a thermodynamic phase transition, from both the asymmetric 
sponge and the vesicles and disks regions, in keeping with 
the continuity between the Higgs and confinement regions. 
However, the appearance of the surfaces nonetheless varies 
drastically. Deep in the Higgs phase, long edges frame large 
surfaces. By contrast, in the 'seaweed' regime, edges are still 
cheap and the tradeoff is between edge entropy and the cost of 
the surfaces attached to the edges, leading to lines joined by 
(lattice dependent) surfaces of close to minimal area. 
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FIG. 2: The line tension/Fredenhagen-Marcu diagnostic is the ratio 
of a half Wilson loop (numerator) to the square root of a Wilson loop. 
Fat points denote matter insertions. On the left we give a cartoon of 
the dominant behavior in the confined phase and on the right the 
dominant behavior in the deconfined phase. 



IV. THE HUSE-LEIBLER HORSESHOE 



properties of surface and edges. Once this has been done, the 
resulting ratio will either vanish as exp(— Xg) in the symmet- 
ric sponge, or approach a constant. Indeed, one does not even 
need to introduce a separate length scale g altogether. One can 
instead decide to cut the Wilson loop in two. For a square loop 
of side length L, as shown in Fig. [2] this yields two rectangular 
loops of size L/2 x L, one long side of which is open and now 
takes over the role of the gap; again, a x needs to be attached 
to each open end. The result is the line tension diagnostic ratio 
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where the curves C l /i an d C define the half- Wilson loop Wj / 2 
and full Wilson loop W respectively. R(L) has the desired 
asymptotic behavior: 



limL->ce R(L) = deconfined phase 
limi^oo^(L) t^O otherwise. 
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Having described the phase diagram of the Z 2 gauge theory 
in terms of our system of surfaces with edges, let us now turn 
to the task of distinguishing the symmetric sponge/deconfined 
phase from all others. Intuitively, we need a quantity which 
teases out whether the perimeter law of the Wilson loop is 
fundamentally due to a vanishing surface tension, or instead 
due the absence of a surface spanning it and whether there 
is an underlying large line tension. Somewhat non-intuitively 
this can be accomplished by simply defining an appropriate 
line tension. 

To this end consider the geometrical object inscribed in 
Fig-Hi It is a Wilson loop broken open, so that its total length 
remains fixed but so that gap of size g opens up. To keep the 
quantity gauge invariant, it is necessary to terminate the open 
strings by a matter term, so that the Huse-Leibler horseshoi 
is defined by the quantity 
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H(L, 




(6) 



where the product is taken over the links Ch (s,s ! ) constituting 
the horseshoe between s and s' . 

The basic point is that, if surfaces are cheap and edges 
are expensive — the situation prevailing in the symmetric 
sponge — the system will cover most of the horseshoe with a 
surface and then close the gap with an edge at a cost in free 
energy of exp(— Xg) due to the tension of the line closing the 
gap (red line in Fig.Q]). By contrast, if surfaces are expensive 
or/and edges are cheap which is the situation elsewhere in the 
phase diagram, no such dependence on g arises. As a practical 
matter, to separate out the dependence on g from that on the 
linear size of the horseshoe, L, we need to study this quantity 
in the limit L ^> g ^> 1 and then check for an effective line 
tension X in the expectation value of the horseshoe. 

An elegant alternative is to remove the L-dependence by di- 
viding the horseshoe by the Wilson loop of same size L, as the 
leading L dependence of the two is set by the thermodynamic 



Before we work further with R, three comments are in or- 
der. First, along the K = °° axis this ratio is simply the spin- 
spin correlation function of the Ising model evaluated between 
the endpoints of the half Wilson loop and it distinguishes the 
two phases as asserted. Second, this diagnostic fails along ex- 
actly one line — when J = 0, but fortunately the Wilson loop 
itself is available exactly in that case. Third, we have checked 
by Monte Carlo and in the small J expansion, that R has the 
claimed behavior. 



V. OPERATOR FORMULATIONS AND THE 
FREDENHAGEN-MARCU ORDER PARAMETER 

The line tension diagnostic introduced in the last section 
can be oriented in an arbitrary direction for the Lorentz invari- 
ant example studied thus far. To obtain an operator formula- 
tion, we must pick a particular orientation. Here we consider 
three and gain some insight into the operation of the diagnos- 
tic. 



A. Equal time formulation 

The simplest option is to orient the variables in the ratio so 
that they all lie in the same time slice of the path integral as 
in Fig. |3(a)| Then the transcription to the Hamiltonian formu- 
lation is straightforward and we obtain from the definition in 
Eq. 0: 



, x Wu 2 {L) 
R(L) = 
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where \G) is the ground state. As written, at T = 0, R(L) 
clearly measures a property of the ground state wavefunction 
but it is straightforward to use it at T > with thermal aver- 
ages replacing ground state ones. Transcribing our previous 
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analysis, R(°°) vanishes if the equilibrium gauge field fluctu- 
ations are weak and the matter is uncondensed but not other- 
wise. This turns out to be the form of the diagnostic which 
can be used most generally for reasons that we explain later. 



B. Fredenhagen-Marcu Order Parameter 

Now let us orient the operators so that they lie in a plane 
containing the imaginary time axis with the endpoints of the 
half Wilson loop at the same time (see Fig. |3(b)| i which we 
take to be zero for convenience. Then we can write the opera- 
tor expression 

W 1/Z =(G\^ S , n°f(-r/2)|G) (10) 

where T = cL up to multiplication by a velocity c, here set to 
1, C ss i is the linear path between s and s', and 



af(-T/2) = e 
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and e is the transfer matrix for one step in the time 
direction — with the implied meaning in the time continuum 
limit. [In the process of passing to the Hamiltonian setting we 
have, as usual, set a z = 1 along the temporal links.] 
Similarly, we can write the denominator as 



W = (G\ U&j(+T/2) n°f(-r/2)|G> 

iec.j lec.j 
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(12) 
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as the candidate two spinon state constructed from the ground 
state. Thus we see that 



R(L) 



(G\ss') 



(14) 



i.e. it is the overlap between the ground state and the normal- 
ized two spinon state. In this form, our diagnostic was first 
discovered by Fredenhagen and Marcu 13,14 and they proved 
that it does show the asymptotic behavior ^ in the Z2 gauge 
theory with matter. In subsequent work they also considered 
the equal time form (O. 

It follows from our previous considerations that 

• in the deconfined phase the candidate two spinon state is 
orthogonal to the ground state in the limit cL = T — >• 00, 
and has finite energy and thus there are free spinons in 
the spectrum 

• in the confined-Higgs phase the candidate two spinon 
state is not orthogonal to the ground state and thus there 
are no free spinons in the spectrum 

In this fashion we see that in this interpretation, R(°°) dis- 
tinguishes between the phases based on the existence of free 
spinons in the spectrum. 

Finally, it is instructive to see how these limits work directly 
in the Hamiltonian formalism. Consider the sequence 



(a) the gauge-noninvariant string YheC 1 °/( — ^/^) acts on 
the ground state leading to a state with a modified con- 
straint at s and s 1 

(b) the system evolves/relaxes for an imaginary time T /2 

(c) the gauge-noninvariant product xfTy acts on the state 
at the end of the last step and leads to a final gauge 
invariant state 

Deep in the deconfined phase (a) and (b) leave an additional 
electric flux string that meanders between the sites s and s' 
and (c) adds the matter needed to create a proper two spinon 
state. Conversely, deep in the confined phase, (b) eliminates 
the string introduced by (a) in favor of fluxless charges at s 
and s' which are then eliminated in (c) leaving the ground state 
behind. 




(a) Equal time diagnostic 




(b) Fredenhagen-Marcu "Order Parameter" 




(c) Spinon derealization diagnostic 



FIG. 3: Different operator formulations of the deconfinement diag- 
nostic ifTJ. The blue curve and terminating points define the locus of 
numerator and the green curve marks the needed completion for the 
denominator. 



C. Spinon derealization 

The two interpretations provided above are sufficient for 
our purposes in the rest of the paper, but some further insight 
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can be gained by considering a third one: Now we orient the 
half Wilson loop so that s and s' are separated by imaginary 
time T as shown in Fig. |3(c)| 

With this choice, the numerator becomes the expectation 
value, 

W l/2 = (G\zl(T/2) J] af(772K(-T/2) ]J °i(-T/2)\G) 

lec 0s , leC 0s 

(15) 

while the denominator is given by the same expression exhib- 
ited previously. 

To understand why the ratio W\n/\W behaves differently 
in the two phases consider the intermediate states produced by 
imaginary time evolution. In the deconfined phase 

• ^1 12 couples to a state with one defect site at s' and a 
delocalized spinon, while 

• W couples to a state with two widely separated defect 
sites. 

Hence the ratio is controlled by intermediate states with dif- 
ferent energies and scales as 

g~ (^defect^^spinon)-^ 

== (16) 

V e — (£defect+£defect ) T 

and thus vanishes at large cL = T. 

By contrast, in the confined phase, both operators couple to 
eigenstates with localized spinons, at s' for W\ h and at both 
s, s' for W, whence the ratio tends to a constant. It follows 
then, that in this version the diagnostic is sensitive to the exis- 
tence of a delocalized spinon state in the deconfined phase. 

VI. FLUCTUATING CONSTRAINTS AND CHOICE OF 
OPERATOR FORMULATION 

In a condensed matter context, a constraint such as the 
one encoded in Eq. [2] is not microscopic (or "fundamental") 
but rather of emergent origin, that is to say it is a legacy of 
higher energy terms in the Hamiltonian. The constraint is thus 
not inviolable — rather, states outside the "gauge invariant" or 
"physical" sector are associated with a large but finite energy 
penalty set by a scale U. To put it differently, in the context of 
gauge theory it is assumed that "unphysical" states are really 
unphysical — they are forced on us by the challenge of han- 
dling the constraint — but in condensed matter setting they are 
certainly physical but are instead high energy states. 

This distinction has an important consequence for the oper- 
ator versions of our diagnostic, e.g. the Fredenhagen-Marcu 
version (fT4-b . As noted in the discussion following its intro- 
duction, it contains a two spinon state creation operator that 
can be factored into two pieces, 

tft and n °l(-T/2) 

which are not individually gauge invariant, although their 
combination is. Now the second piece is sensitive not only 



to matrix elements of H between gauge invariant states but 
also to matrix elements of H involving gauge non-invariant 
states. Indeed our verbal explanation of the creation of the 
two spinon state referred specifically to imaginary time evo- 
lution in the gauge non-invariant sector where the constraint 
was violated at two sites. For the specific choice of Hamilto- 
nian (Q]) all is well — the detour into the gauge non-invariant 
sector introduces no physics inconsistent with the proper op- 
eration of the diagnostic. Of course this had to be the case as 
we got this form working from a manifestly Lorentz invari- 
ant Euclidean path integral. We now show that this does not 
have to be the case by explicitly exhibiting an example with a 
soft constraint where the Fredenhagen-Marcu version of R(L) 
fails to diagnose deconfinement. The example is somewhat 
contrived but it is intended only to function as a cautionary 
counter-example . 

Consider modifying our starting Z2 gauge theory Hamilto- 
nian (fTJ to Hf. 

H f= H + T E(G.v - l) 2 - (G s 1)t* (17) 

and set i,r->0 for now. Now we do not implement the con- 
straint © on the Hilbert space, instead a suitably large value 
of U does this in the low energy sector. With that ordering of 
energy scales, the ground state at small X <C U 

1 c ) = (n \ (g* - 1)) 1 k = id® 1 k = id (is) 

has precisely the same form as that of the toric code. However, 
there is now a change of meaning, as microscopically both 
the link and site variables are physical while in the toric code 
the site variables can be gauged away. Yet, at low energies 
in our model the same reduction takes place and so it has as 
many degrees of freedom and upon our choices of parameters 
exhibits the same, but now emergent, Z2 order. 

Now consider the excitations of our model. The gauge in- 
variant excitations preserve G s = 1 at all sites and are visons 
with an energy 2K and spinons with an energy 2Fm- Note 
that their energies are independent of X; indeed, all gauge in- 
variant states have (by construction) energies independent of 
X. The gauge non-invariant excitations involve violating the 
constraint by setting G s = — 1, and also come in two species: 
spinonic defects with energy U + 2Fm — 2X obtained by set- 
ting x x = — 1 at a site and gauge defects with energy U + 2X 
obtained by acting with a string Yl^i extending from a given 
site to infinity. Observe that the energies of the two defect 
states cross when U ^ X > Ym /2 even though the spectrum 
of the entire set of lower energy gauge invariant states is un- 
changed. [For the Hamiltonian ([T]) the spinonic and gauge 
defects have energy 2Tm and respectively. While there is 
no explicit energy penalty for gauge non-invariant states, the 
separate requirement that the ground state obey (f2]i eliminates 
the degeneracy that would be present otherwise.] 

In our previous discussion we argued that the action of 
ll/ec /°/( — T /2) was to create a string stretching between 
sites s and s' . This works as long as the state with two gauge 
defects is the ground state in the sector with two defect sites. 
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Once the defects cross and the two spinonic defects become 
the ground state, this is no longer true. As a technical matter, 
the crossing alone is not enough, we also need to turn / on 
to a small value so that imaginary time evolution connects the 
two states. Once this is done, R(°°) 7^ even though we are 
clearly in the deconfined phase and the diagnostic fails in this 
form. 

However the diagnostic still works in its equal time form 
so the lesson is that in general condensed matter settings the 
formulations presented in the previous section are not fully 
equivalent: whereas the equal-time formulation (Sec. IV A| ) is 
robust, the other pair of formulations (Sec. IVBIVCI ) need not 
be. 



VII. FINITE TEMPERATURE TOPOLOGICAL ORDER IN 

d>3 

We now return to our main development and consider the 
use of our diagnostic to establish the finite temperature phase 
structure of the Z2 gauge theory where contrasting answers 
can be found even in the recent literatur o 15 ' 18 ' 21 . We begin by 
considering the T > physics of the toric code, J = T = in 
Eq.Q] It can be shown (see Appendix [Ab that at this point in 
d + 1 space-time dimensions the quantum partition function at 
all temperatures is proportional, up to a simple multiplicative 
analytic factor, to that of the pure Z2 gauge theory without 
matter defined on a d dimensional Euclidean lattice. 

From what is known about the latter theory, one can im- 
mediately read off that topological order will persist to fi- 
nite temperature for all d > 3, but not below. Further we see 
that, despite the presence of matter, the transition between the 
low temperature deconfined phase the the high temperature 
phase can be diagnosed by the Wilson loop which goes from 
a perimeter law to an area law between the two phases. 

This physics has to do with the gauge excitations. For 
d = 2, finite energy point-like visons appear at an exponen- 
tially small but non-zero density at any T > destroying the 
deconfined phase. As the Wilson loop measures the averaged 
parity of the enclosed vorticity, it acquires an area law as a 
result. 

By contrast, this mechanism is not operative in higher 
dimension because the gauge/magnetic excitations now ap- 
pear in a higher-dimensional version — loops in d = 3 — whose 
characteristic size, i; m , vanishes as T — > 0. As the Wilson loop 
now counts the parity of flux loops which link it, it is sensitive 
only to those within £,„, of its boundary, whence the perime- 
ter law and the deconfined phase survive at sufficiently small 
temperatures. 

At higher temperatures there is a phase transition, e.g in 
d = 3 where the vortex loops unbind. The basic physics of 
vortex loops discussed above indicates that the finite temper- 
ature deconfined phase should be stable beyond the Kitaev 
point in d > 3 as argued i n 15 ' 18 . However, once matter be- 
comes dynamical the Wilson loop exhibits a perimeter law at 
any T and other standard diagnostics fail as well. 

Fortunately, our diagnostic continues to work at T > 0. At 
small temperatures one can work perturbatively at small J, 



when the dynamical matter is heavy. As at T = 0, the basic 
result is that in this limit the behavior of R(L) follows from 
that of W(L) in the absence of dynamical matter. As the latter 
exhibits a transition between perimeter and area dependence, 
R(°°) transitions between vanishing and developing an expec- 
tation value of 0(L°). In somewhat more detail, one can ex- 
amine the Euclidean path integral, now finite in the (large) 
imaginary time direction at low temperatures much as in Sec- 
tion |nil As before, a large line-tension between any pair of 
charges survives, and it continues to be preferable to mini- 
mize free edges at the expense of adding the (cheap) surface, 
leading to a vanishing of R(°°). 

At high temperatures, T — > °°, it is more convenient to 
work directly in the Hamiltonian formulation and the high- 
temperature expansion. In this expansion, one obtains poly- 
nomials in powers of |3 = l/T upon expanding exp(— $H). 
Non-vanishing terms in the expectation values appearing in 
Eq.[7]are only obtained if each unpaired o : ,X : is matched by 
the polynomials. The cheapest (in powers of p) way of doing 
this is to retrace the perimeter of W and W1/2. This leads to 
the same leading power of p in the expectation value of each, 
which thence cancel in the expression for R(°°). Consequently 
^?(°°) is non-zero at sufficiently high temperatures and we that 
there is a phase transition between the high and low tempera- 
ture phases. 

VIII. BEYOND Z 2 
A. Gauge fields with fundamental matter 

Thus far our discussion has focused on the simplest case, 
that of standard Z2 gauge theory with Z2 matter where there 
are two phases to distinguish — the deconfined phase and the 
confined-Higgs phase. Our deconfinement diagnostic general- 
izes straightforwardly to lattice gauge theories based on other 
gauge groups where there are similarly only two phases to 
distinguish, to Z„, 1/(1) and to Yang-Mills theories based on 
non-abelian groups. Indeed, the last named case was the pri- 
mary spur to Fredenhagen and Marcu's work. In all of these 
cases R has the same structure with appropriate replacement 
of the matter and gauge fields. 

For example, a compact U(l) gauge theory with charge q = 
1 matter is defined by a space time action analogous to (0]i 

-S = K^Y\Uj {l) +J £ T^v+h.c, (19) 

P ledp l;s/edl 

where the Ui = exp(/A/) = £/*, are now 1/(1) valued gauge 
fields living on the links defined in the standard orientation 
from one sublattice to the other, and the x s = exp(/0 s ) are the 
corresponding matter fields residing on the sites. Unlike the 
Z2 theory, we now need to pay attention to the orientation of 
various link variables. The links in the second term are tra- 
versed in the standard orientation while in the first term the 
loops are taken anticlockwise with /(/) = ±1 if the link is tra- 
versed in, or opposite to, its standard orientation respectively. 

At zero temperature in spatial dimension d>3 this theory 
exhibits both a confined and a deconfined phased see Fig. [4] 
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Absent matter, 7 = 0, the oriented Wilson loop 

wi[c] = ([\u m ) 

lec 



(20) 



is capable of distinguishing between these two phases, ex- 
hibiting an area law in the former, and a perimeter law in 
the latter, with logarithmic corrections signaling the Coulomb 



interaction 



20 





Higgs 


1^(00) *0 




Wi - e _L 


deconfined 


confined / 


R^oo) = 




Wi - e" L 



K 



FIG. 4: Sketch of the phase diagram of the f/(l) gauge theory with 
charge one matter fields. The behaviour of the line tension diagnostic 
R\ and Wilson loop W\ is indicated for each phase. 

This phase structure persists when matter with charge q = 1 
is turned on. However, as in the case of the Z2 theory above, 
the Wilson loop fails as a diagnostic in this setting. By anal- 
ogy to Eq. [7] we construct the appropriately generalized diag- 
nostic ratio, 



«i(L) 



W l/2 (L) _ (^(Uiec 1/2 Ul(,)^') 



(21) 



As before the phases are distinguished by the asymptotic be- 
havior of R\(L), 



limi^oo R\(L) =0 deconfined phase 
hm^oo^i(L) 7^0 otherwise. 



(22) 



We now turn to the application of the diagnostic to cases 
with more than one phase and to problems phrased in more 
standard condensed matter settings. 



B. U(l) gauge theory with charge q matter 

To obtain a richer phase structure we consider a £7(1) gauge 
theory coupled to matter fields that carry one or more higher 
(non-fundamental) charges q > 2. In this setting it is useful to 
consider the set of charge-^ operators 



lec 



(I)' 



R q (L) = 



Hlec%)> 



(23) 



whose structure indicates that q electric flux quanta emanate 
from a given charge. 





Higgs, Z 2 deconfined 


confined 






y r 2 (00)^0 




\ W 1>2 - ei>- — " 


R 2 (°°) *Q 




Wi - e _L2 


S R 2 (00) = 


W 2 ~e- L / 


Wi >2 ~ e" L 




deconfined 



K 



FIG. 5: Sketch of the phase diagram of the 1/(1) gauge theory with 
charge two matter fields. Charge-g line tension diagnostic R q and 
Wilson loop W q are defined in Eq. 1231 



As our first example we consider the well known phase di- 
agram of a single matter field of charge q = 2. This contains, 
as shown in Fig. El m addition to the deconfined phase of the 
f/(l) field, a Higgs phase which corresponds to the decon- 
fined phase of the appropriate Z2 gauge theory and a distinct 
completely confined phase- . Due to the lack of q = 1 matter, 
we can make some progress using Wilson loops alone. 

Of the three phases, two (1/(1) and Z2 deconfined) exhibit 
a perimeter law in W\ while one (1/(1) confined) exhibits an 
area law. The key point is that the q = 2 dynamical matter can- 
not break up the confining string that stretches between two 
test q = 1 sources in the confined phase. By contrast, W2 ex- 
hibits a perimeter law everywhere and distinguishing the two 
deconfined phases now requires that we evaluate the diagnos- 
tic R2- As this tests for the presence of doubly charged states 
in the spectrum, it vanishes only in the fully deconfined phase. 

In terms of the nature of gauge-field fluctuations, these dis- 
tinct behaviours can be seen most intuitively by gauge-fixing 
(j) s = 0. Then, deep in the deconfined phase, an appropriate 
projection of the state <8)/|A/ = 0) is a ground state. In contrast, 
the gauge field fluctuates strongly in the confined phase, and 
the A/ range over the full interval [0,27t).In the Higgs fluctua- 
tions are restricted: here, A/ = or A/ = K satisfy the second 
term of (TT9T>. By construction, R2 is sensitive to the former, but 
not the latter, fluctuations, and it can thus diagnose the Higgs- 
confinement transition. This state of affairs is summarized in 
TableU 



Phase 



Wi(L) R 2 (L) 



Deconfined Perimeter 
Higgs Perimeter L 
Confined Area L° 



TABLE I: Diagnostics for U(l) gauge theory with only q = 2 matter 

As a second example we consider the more generic and 
complex situation where charges with both q = 1 and q = 2 
are present simultaneously. The number of parameters in the 
Hamiltonian makes for a large phase diagram, and we restrict 
ourselves to the instructive case where the q = 1 matter is very 
heavy, i.e. very sparse. This then has the same physics of our 
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FIG. 6: Sketch of the phase diagram of the 1/(1) gauge theory with 
charge two as well as heavy charge one matter fields. 



last example, and thus the same phase diagram as Fig. [5] The 
new feature is that now W\ obeys a perimeter law in all three 
phases so that we may no longer use it to test for the decon- 
finement of charge- 1 objects. The solution to this is to evalu- 
ate R\ instead. The diagnostic behavior of the combination of 
Ri and R2 is summarized in Table iHl 



IX. CONCLUSION 

We have substantially met the challenge posed at the start of 
this paper, that of finding a diagnostic for deconfinement that 
has the form of a ground state expectation value, i.e. which is 
a property of a single ground state of the system. It is also nice 
to find that our reasoning, which comes from thinking about 
line and surface tensions in the statistical mechanics of sur- 
faces and edges, leads to an object discovered by lattice gauge 
theorists thinking about creating states with (non) deconfined 
quarks. 

The next steps are to apply this formulation to problems 
other than those addressed or reviewed in this paper and two 
come to mind. The first is using it to detect spin liquids in 
systems with Heisenberg symmetry where the gauge field in- 
volves valence bonds and there is a set of orthogonality issues 
to contend with. The second is generalize the construction to 
non-abelian topological phases captured by the models con- 
structed by Levin and Wen^E. 

X. ACKNOWLEDGEMENTS 



Phase Rj(L) R 2 {L) 



Deconfined 





Higgs L° 





Confined L° 





TABLE H: Diagnostics for 1/(1) gauge theory with q = 1 and q = 2 
matter 

We note that this combination manages to distinguish all 
three phases, whereas the fail to do so. 

These deconfined phases persist to finite temperatures in di- 
mension d > 4. The known behavior in lower dimensions is 
more complicated. The U(l) deconfined phase is absent in 
d = 3 at T > and in d = 2 even at T = while the Z2 decon- 
fined phase is present in both these cases. This complexity is 
correctly captured by our doublet of diagnostic ratios. 



C. Gauge theories of quantum magnets 

The formalism developed above can be directly applied to 
two classes of Hamiltonians which arise naturally in discus- 
sions of quantum magnetism. The first are the "slave particle" 
or Schwingerboson/fermionHamiltonians22r— . Here the spin 
Hamiltonian of a quantum magnet is reformulated as a matter 
gauge system with spinons coupled to a gauge field. Here our 
construction of R(L) carries over pretty much directly. 

The second class are quantum dimer models 25,2 ^, now with 
a Hilbert space expanded to include gapped, spin-carrying 
monomers to represent spinons 27 . In this setting, the micro- 
scopic gauge structure arises from the hard core dimer con- 
straint and the Wilson loop is defined as a dimer rearrange- 
ment on a loop that respects the hard core constraint 2 ^. Ac- 
cordingly we can again define our diagnostic. 
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Appendix A: Kitaev's toric code 

The model describing Kitaev's toric code is obtained from 
the Ising gauge Hamiltonian, Eq. Q] straightforwardly. First, 
set r = / = 0. Secondly, replace the matter term proportional 
to by utilizing the constraint, Eq. [2] 

G s =l=^x:= ]~[ of (Al) 
ksedl 

to obtain 

-h tc = n °-/+ r «i n °? ■ < A2 > 

p ledp s t.sedl 

which now acts in the subspace spanned solely by the link 
degrees of freedom. 

One can see that the model is special by noting that all 
terms commute with all others. This allows the full spec- 
trum to be readily obtained and one sees that the model ex- 
hibits topological order. This is most easily done before 
gauge fixing, (Eq. lAll ). in the basis of eigenstates of {of} and 
{tJ}. In this formulation, sites where tJ = — 1 are occupied 
by "charges" and links where of = — 1 carry "electric flux": 
G s = 1 is the Z2 version of Gauss's law. The dual, magnetic, 
flux passing through loops C is measured by the Wilson loops 

W[C}=niectf- 

The ground states of Htc lack charges, tJ = 1 , V*, and are 
characterized by flux expulsion, Yliedp®] = 1> Vf- 

The complete lack of gauge field fluctuations, diagnosed by 
the Wilson loop W = 1, signals a deconfined phase. Consis- 
tent with this ground state signature, the excited states contain 
non-interacting charges, free to sit at arbitrary locations. 
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1. Mapping to classical Z2 gauge theory and T > phase 
transition 

We now turn to a computation of the partition function of 
Kitaev's model at T > 0. As gauge and matter sectors of the 
theory only interact via the constraint, G s = 1, one can use the 
basis states |{c/}) ® |{x s }) where o^,|{o/}) = 0/'|{o/}) and 
TylM) =VK X *}) to obtain 

Z = Z L Zs (A3) 

where 

Z L = £ e P^I;.n; e3; ,o, (A4) 

0;=±1 

is the partition function of the classical or discretized Eu- 
clidean lattice Z? gauge theory in d dimensions and 

Zs= E Sn^.ieP^ 1 ' (A5) 

T s =±l 

is the partition function of independent spins in a magnetic 
field r,„ up to a single, global, constraint that the number of 
flipped spins (in gauge language: charges) be even. All the 
singularities in Z thus are contained entirely in Zl. 

The finite temperature phase structure of Zl is well known. 
In d > 3 it exhibits two phases. The deconfined phase, at large 
$K, is characterized by a perimeter law for the Wilson loop at 
large loop sizes. The confined phase, at small $K, exhibits an 
area law. By contrast, in d — 2, Zl exhibits only a confined 
phase, but with a confinement scale (area law coefficient) Aq 
that diverges exponentially as $K —> °°. 

All of this ties in with our previous results on Ising gauge 
theories: (i) Kitaev's model exhibits a T > continuation of 
the T = topological phase in d > 3 but not in d = 2. (ii) 
the presence of the topological phase can be detected by the 
(spatial) Wilson loop exhibiting a perimeter law. 



local perturbation V to the toric code, where locality now im- 
plies that the range of variables coupled is strictly bounded, 
say by a distance b. For example, we may reintroduce the two 
terms in Eq. Q] which we had set to zero on our way from the 
Z2 gauge theory to the toric code, although the argument to 
follow is completely general. 

The easy part is to notice that the topologically distinct, 
degenerate ground states of Hq will not mix up to an order 
L/b ~ 0(L) since they differ by the creation of a vortex loop 
(or appropriate dimensional analog) that has to be created and 
stretched across the system in as many applications of the per- 
turbation. The less obvious part is that the states remain de- 
generate in energy up to that order. 

To see this, observe that the energy shifts &E for two ground 
states, |\|/o) = \G) (Eq.[[8])and =4>- r |\(/o) differing by the 
flux through a non-contractible loop, C nc , are sums of terms of 
the form 

8E = (yo\VPV...\y ) (A6) 
8£i = (\\ro\& c VPV...<f> x \\\fo) • (A7) 

Here V is the local perturbation, P = 
(l-L l¥i)(¥/|)i^bfo (1-L l¥i)(¥/|) projects onto the 
excited state manifold and <£>* = fltec °? inserts a flux by 
operating on all links crossed by the non-contractible loop, 

Now, we claim that 8Eq = 8E1 to o(L) in perturbation the- 
ory. To see this, note that we can permute one <t>- T through the 
operator products of Vs and Ps to annihilate the other. This is 
possible a) because 4^ commutes with both numerator and 
denominator of P and b) because any topologically equiv- 
alent deformation of the non-contractible loop along which 
the <T v s act produces the same state. Thus at an order below 
L/b ~ 0(L) in perturbation theory, there is always a choice of 
loop which avoids the locations of all the perturbing terms. 



2. Stability of topological order under perturbations of the 
toric code 

It is often convenient to do calculations and point-of- 
principle demonstrations for topological phases at a "Kitaev 
point", where terms generically present in the Hamiltonian are 
set to vanish by hand, leaving behind a soluble model. This 
has some obvious drawbacks — for instance, a "zero-law" for 
the Wilson loop, W = 1 = e~ L is obviously somewhat non- 
generic, like a vanishing rather than finite correlation length 
in a disordered system, and any result thus obtained needs to 
be supplemented at least by a discussion of its fate under per- 
turbations. 

We now show that the deconfined phase is characterized 
by a persistence of the 2 A '" 1 -fold topological degeneracy that 
signals topological order, where N nc is the number of inde- 
pendent non-contractible loops. More precisely, it exhibits a 
cluster of states with splitting of 0(e~ L ) for a system of lin- 
ear dimension L which are separated from all other states by 
a gap of 0(L°). To this end, we consider adding an arbitrary 
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